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1 Introduction 

In one dimension the so-called Riccati equation takes the form 

This non-linear equation is also called Miura transformation and gives a 
relation between the Korteweg-de Vries equation and the modified Korteweg- 
de Vries equation. It is related to the one-dimensional Schrodinger equation 

-1^—0 (2) 

with V a function called the potential, in the following way. Given a particular 
solution for ([2]), one can construct the corresponding solution for ([T]) simply 
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by taking the logarithmic derivative y = — — and vice versa. 

u ax 

Moreover, there is a second relation between the one- dimensional Schrodinger 
equation and the Riccati equation, namely the one-dimensional Schrodinger 
operator can be factorized in the form 

dx"^ ^dx ^dx ^) ^ 

if and only if y is a solution of ([1]). 

An analogous factorization of the three-dimensional Schrodinger operator 
is obtained in |2] (see also |5]) using quaternionic analysis, while some re- 
presentation formulae for the solution of the three-dimensional Schrodinger 
equation are contained in pQ. 

The Riccati equation has some peculiar properties among which two theo- 
rems of Euler. The first of these states that given a particular solution of 
the Riccati equation, the general solution can be found in two integrations; 
while the second one states that given two particular solutions, the general 
solution can be found in one integration. In [H] a quaternionic generalization 
of the Riccati equation is established and versions of the above-mentioned 
Euler theorems corresponding to this generalization are proved. 

Closely connected to the Schrodinger operator is the Darboux transfor- 
mation. In general it can be defined as follows. Consider the eigenvalue 
equation 

Lf = \'f (4) 

with A a complex number. Assume that the operator L admits the factori- 
zation L = AB with B a linear operator. Then naturally dl]) yields 

Lifi = \'fi , 

where Li = BA and /i = Bf. Thus applying B to eigenfunctions of L yields 
eigenfunctions of Li. 

Although this scheme is general, the Darboux transformation has been stu- 
died mostly for the Schrodinger operator. In view of the factorization ([3]), 
the Schrodinger equation 

dx'^ ) 
has the following Darboux image 

Lifi = -(^-y) (^ + y)fi = >^'h 

\dx J \dx J 
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where /i = ( y ] f- The operator Li has again the form of a Schrodinger 



operator 



dx 



^ dx^ ^ 



1 dy 2 

where Vi = — y . 

dx 

This paper is devoted to a factorization of the higher dimensional Schro- 
dinger operator in the framework of ChfFord analysis, a direct and elegant 
generalization to higher dimension of the theory of holomorphic functions in 
the complex plane. This factorization is based on one of the most fundamen- 
tal features in Clifford analysis, viz. the factorization of the Laplace operator 
A„ by the so-called Dirac operator ■ — A„ = d"^. It is precisely this Dirac 
operator which underlies the notion of monogenicity of a function, a no- 
tion which is the multi-dimensional counterpart to that of holomorphicity in 
the complex plane. 

The outline of the paper is as follows. For the reader who is not familiar 
with Clifford analysis, we recall some of its basics in Section |2J In Section Owe 
establish a factorization of the higher dimensional Schrodinger operator using 
Clifford analysis. As in the one- dimensional case, this Clifford factorization 
of the Schrodinger operator yields a generalization of the classical Riccati 
equation to the Clifford analysis setting. Moreover we show that for various 
reasons this so-called Clifford Riccati equation is a good generalization of 
the classical one. Some possibilities for obtaining particular vector-valued 
solutions of the Clifford Riccati equation are discussed in Section HI Fur- 
thermore, in Section [5] generalizations of Euler's theorems for the Clifford 
Riccati equation are established. Next, a generalized Schrodinger operator 
in Clifford analysis is studied (Section |6]). In a final section a decomposition 
of the kernel of this generalized Schrodinger operator is obtained. 



2 Some basic notions of Clifford analysis 

Clifford analysis (see e.g. [3] and offers a function theory which is a 
higher dimensional analogue of the theory of the holomorphic functions of 
one complex variable. 

The functions considered are defined in the Euclidean space M" {n > 
1) and take their values in the Clifford algebra Mo,n or its complexification 
Cn = Mo,n ® C If (ci, . . . , e„) is an orthonormal basis of M", then a basis 
for the Clifford algebra Ro,n is given by (e^ : A C {I, . . . , n}) where 60 = 1 
is the identity element. The non-commutative multiplication in the Clifford 
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algebra is governed by the rules: 

CjCik + cikCj = -2 5j^k , j,k = 1,2, ... ,n . 

Conjugation is defined as the anti-involution for which 

Cj — —Cj , j — 1,2, . . . ,n 

with the additional rule i — —i in the case of . 
For k — 0,1, . . . ,n fixed, we call 

mJ^ = |a e Ro,n ■■ a=J2 ^^^^ ; «a e m| 

\A\=k ^ 

the subspace of /;;- vectors, i.e. the space spanned by the products of k different 
basis vectors. The 0-vectors and 1-vectors are simply called scalars and 
vectors respectively, while the vector cat = 6162 ... e„ is called the pseudo- 
scalar. 

The Euclidean space is embedded in the Clifford algebras ]Ro,n ^-nd C„ 
by identifying {xi, . . . , x„) with the vector variable x given by 



2=E 



GjXj 



The product of two vectors splits up into a scalar part and a 2-vector or a 
so-called bivector part: 

xy = x • y + xAy , 

where 

n 

X • y^-{x,y)^-Y^ Xjyj 

and 

n n 
i=l j=i+l 

Note that the square of a vector variable x is scalar-valued and is the norm 
squared up to a minus sign: 

^ — ~ ■ 

The eUiptic vector differential operator of the first order 



called Dirac operator, splits the Laplace operator in M.^ : 

For a differentiable scalar-valued function and a differentiable Clifford 
algebra- valued function /, we have 

9.(0/) = / + . (5) 

Furthermore, this Leibnitz rule admits the following generalization in the 
case of a A;-vector valued function Gk = J2\A\=k^k,A (see for example [7j) 

n 

d,{Gkf) = d^{Gk) f + 2j2 [ejGk]k-i 9,,(/) + (-l)'G, d,if) (6) 

i=i 

where [a]k denotes the projection of a G Mo.n on Mq^^ . 

3 Clifford factorization of the Schrodinger ope- 
rator and the CHfford Riccati equation 

In the sequel, denotes the operator of multiplication with / from the 
right, i.e. M^^g = gf . 

The following results generalize to higher dimension the relations mentioned 
in the introduction between the one-dimensional Schrodinger equation and 
the Riccati equation. 

Proposition 3.1 The Schrodinger operator may be factorized into 

(-A„ - vl)<p = [d^ + Mf) [d^ -Mf)(t> , 
where I is the identity operator and a scalar-valued function, if and only if 

dj + f = v . (7) 
Proof. Applying the Leibnitz rule yields 

{d, + Mf){d,-Mf)(P = {d, + Mf){d,<P-<Pf) 

= -A„0-4(0)/-0a^(/) + a^(0)/-0f 

= -A„0-(4/ + /2)0 . 

Hence 

(-A„ - t;/)0 = {d^ + Mf) {d^ - Mf) 

is equivalent with 

dj + f^ = v . □ 
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Proposition 3.2 The scalar-valued function (j) is a solution of the Schrddinger 
equation 

An(j) + V(f) = 

if and only if f = -7^ is a solution of 

dj + f = v . 



Proof. Assume that is a solution of the Schrodinger equation. For / 
we then have 



^f^f2 _ (^.0)^ A„0 , (40)^ 
CxJ + J - ,2 ^ + ^2 



= V . 

Conversely, suppose that / satisfies d^f + f^ = v and that there exists a 



function such that / = For this function (j) we then obtain 

(f) 

A„0 ^ (40)^ 

or equivalently 

A„0 + -y^ = . □ 

As Propositions 3.1 and 3.2 are the counterparts to the two relations men- 
tioned in the introduction between the Riccati equation ([1]) and the one- 
dimensional Schrodinger equation ([2]), equation ([7]) can be considered as a 
good generalization of ([T]) to the Clifford analysis setting. Moreover, for a 
one-dimensional solution / = f{xk)ek equation ([7]) reduces to the classical 
one- dimensional Riccati equation 

d.J + f = -v . 

Hence we call (I7j) the Clifford Riccati equation. 

Remark 3.1 In the special case of a vector-valued solution, the Clifford Ri- 
ccati equation is equivalent with a scalar elliptic partial differential equation. 
Indeed, for a vector-valued solution f equation ^ only consists of the follo- 
wing scalar and bivector part: 

-{d.J)-{lJ) = V , (8) 
4A/ = . 
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The bivector part 

dxjj -d^Ji = forz<j 

implies that for a simply connected domain Q in M", there exists a scalar- 
valued function (p such that 

/ = V0 = dr,, {(p) ei + 9^.2 (0) 62 + . . . + <9^„ {(f)) Cn . 
Inserting the above in the scalar part ^ yields 

An(j) + {V(j), V0) = -V . 

Remark 3.2 From Proposition 3.2 it is clear that vector-valued solutions of 
the homogeneous Clifford Riccati equation 

dj + f = 

d (p 

take the form f = with (f> a harmonic function, i.e. (j) E ker(A„). 

4 Particular vector- valued solutions of the Cli- 
fford Riccati equation 

In this section we discuss some possibilities for obtaining particular vector- 
valued solutions of the Clifford Riccati equation ([7]). 

First we consider the special case in which the potential v takes the form 

V{X) = Vi{xi) + V2{X2) + . . . + Vn{Xn) , 

viz. the variables can be separated. Assuming that 

/ = /i(a;i)ei + f2{x2)e2 + • • • + /„(a;„)e„ , 

the Clifford Riccati equation ([7]) reduces to the system of ordinary one- 
dimensional Riccati equations 

dxjk + fl = -Vk , k = l,2,...,n . 

Hence in this special particular vector-valued solution of ([7]) can be 

found if and only if each of the above one-dimensional equations can be 
solved. 

By means of the existence of a large class of vector-valued solutions of 
the homogeneous Clifford Riccati equation (see Remark 3.2), we are able 
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to reduce the Clifford Riccati equation ([7]) to a scalar differential equation. 
Indeed, suppose that / and g are two vector-valued solutions of the homo- 
geneous Clifford Riccati equation, i.e. 

where 0i and 02 are harmonic. Then the sum / + (7 is a solution of ([7]) if and 
only if we have 

or equivalently 

In particular, choosing 0i = 02 = 0, equation ([9]) becomes 



Hence if G ker(A„) is a solution of the above equation, then 2^— is a 



solution of the Clifford Riccati equation. 



5 Generalizations of Euler's theorems for the 
Chfford Riccati equation 

In this section we establish generalizations of Euler's theorems mentioned in 
the introduction for the Clifford Riccati equation. 

Proposition 5.1 Let h be a particular vector-valued solution of the Clifford 
Riccati equation Then the vector-valued function 

L = 9 + h (10) 

,s also a soluUon of .fun, only ,f , = ^ /„. son., saiur-valuei soluUon 

of the equation 

A„0 + 2(a^0,/i) = . 
Proof. Substituting ffTOl) in ([7]) gives 

d.g_-2{g_,h)+g^ = ^ . (11) 
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As in Remark 3.1, the bivector part of flTT]) implies that for a simply connected 
domain Q in R", there exists a scalar- valued function (pi such that 

g = d^cpi . 

Putting (p = exp{(f)i), we then have that 

dx(p 
(j) 

and equation (ITT]) written in terms of (p becomes 

{d^cPf A„0 2 {d^cPf 

T'^^^^'^'-^^^^^ 

or equivalently 

A„0 + 2(40,/i) = . □ 



Proposition 5.2 Let g = d^cpi, h = dx<p2 be two particular vector-valued 
solutions of the Clifford Riccati equation ([^. Then 

a q — h 

f = ^^ 
— a — 1 

is also a vector-valued solution of where a = K exp{(pi — (p2), K E C 

Proof. A straightforward computation yields 

d^{a){ag-h) dj^a)g + ad^g - d^h {ag - Kf 
dj + f = -, =T^ — + - - ■ - 



—a{g — h)"^ + {ag — h)'^ + a{a — l)dx g — {a — l)dx h 

a{a - 1) (a^ g + f)-{a- 1) {d, h + if) 
(a - 1)2 



where we have used the assumption that g and h are particular solutions of 
the Clifford Riccati equation ([7]). □ 

Remark 5.1 It should he noted that the new solution obtained in Proposition 
5.2 is not necessarily the general solution for vector-valued functions. Indeed, 
let n > 2 and consider the Clifford Riccati equation with v = —1 

dj + f = -l . 
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Two vector-valued solutions of the above equation are g = ei = d^Xi and 
h = 62 = dxX2- Therefore, by Proposition 5.2 we get the class of solutions 

eiK exp{xi - X2) - 62 ^ ^ 

1 ^ ^ , A G (L . 

K exp(xi — X2) — 1 

However the solutions 6j (j > 2) are not included in the above expression. 

6 A Generalized Schrodinger operator 

Let 

d[ := + Mf and d^^ := d^- . 
In this section we consider the equation 

d[d-fg = X'g (12) 

where g is a. Chfford algebra-valued function and A 7^ a complex number. 
We assume that / is a solution of the Clifford Riccati equation ([7]). 
From Proposition 3.1 it follows at once that for a scalar- valued function (p 
equation (fT2|l reduces to the Schrodinger equation 

(_A„ - = A^c/, . 

Hence equation f[T^ will be referred to as the generalized Schrodinger equa- 
tion. 

Naturally the Darboux scheme described in the introduction works for 
the generalized Schrodinger operator L = dld~^ . 

Proposition 6.1 If g is a solution of / fl^) then h = g is a solution of 
the equation 

where Li = dl . 

By means of ([6]) we obtain for a /c-vector valued function Gk '■ 

= {d, + Mf){d,Gk-Gkf) 

n 

= -A^Gk - dJ,Gu)f - 2 [e,Gk]k^i <9,,(/) - {-ifGk dj,f) 
+ dJ,Gk)f-Gkf^ 

n 

= -/\nGk + Gk{{-lf^' - f^) -2Y,[efiu]k-id,J (13) 
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and similarly 

n 

d-fd[Gk = -AnGk - Gk ((-1)'+' dj + f)+2j2 [e,Gk]k-i . (14) 

i=i 

Remark 6.1 It is easy to check that for a scalar-valued function (f> we have 

d^fd[<P=-An<P + <PidJ-f') . (15) 

Proposition 6.1 and equation f[T^ lead to the following result. 

Proposition 6.2 Let f be a vector-valued solution of the Clifford Riccati 
equation and (j) a scalar-valued solution of the Schrddinger equation 

{-An-vI)(f> = X^<j) . (16) 

Then the vector-valued function 

n 

h = hjCj = dx~4> 

is a solution of the equation 

n 

h{-An -v)-2j2hj d.,^l = W . 

Proof. As (fT6|) is equivalent with 

dx- (p = , 

Proposition 6.1 implies that the function /i = is a solution of 

dx~ h = X^h . 
In view of fll4p the above equation is equivalent with 

n 

-Anh - h{dxf_ + /') - 2 5^ h, dx^l = W 
which gives the desired result. □ 

Similarly we obtain from Proposition 6.1, equations f|T3|) . f|T4|) and f|T5l) ano- 
ther result. 
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Proposition 6.3 Let f be a vector-valued solution of the equation 

u = dj-f. 

Furthermore assume that g = Yl]=i 9j^j i^ ^ solution of 

n 

g_{-A^ + u) + 2Y,93d.,l=\\ • (17) 
i=i 

Then the scalar-valued function (p and the bivector valued function H2 con- 
structed in the form 

<P = [dx- g]o , H2 = [dx- g\2 

satisfy the equation 

n 

Proof. Equation (|T7I) can be rewritten as 

dx- g = >?g_ ■ 

Applying the Darboux scheme yields that the function h = dx~ g , which 
consists only of a scalar and a bivector part, satisfies the equation 

dx- dj h = . 

Denoting by and H2 the scalar, respectively bivector, part of h and using 
expressions f|T4|) and f|T5|) leads to the desired equation. □ 

Remark 6.2 Propositions 6.2 and 6.3 are special cases of the following more 
general result. 

Proposition 6.4 Let f be a vector-valued solution of the equation 
Furthermore assume that Gk is a k-vector valued solution of 

n 

Gu{-An + w) -2Y,[(^jGu]k-idx,l = X'Gk . 
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Then the {k — l)-vector valued function H^-i and the {k + l)-vector valued 
function H^^i, constructed in the form 

Hk-i = [dx~Gk]k-i , Hk+i = [dx~Gk]k+i 
satisfy the equation 

(Hk^i + Hk+i)i-An + w) + 2j2 ([ejHk-i]k-2 + [ejHk+i]i}j d^^l 

j=i V / 

7 Decomposition of the kernel of the genera- 
lized Schrodinger operators 

We start this section with another representation of the equations f|T2|) and 

d-fd[h = X^h , Oy^XeC . (18) 

Tl 

Proposition 7.1 Let — be odd and f a vector-valued function. Then 
(i) equation can be rewritten in the form {A^ — \^)g = 0, where 

(a) equation / fi^) can be rewritten in the form [B"^ — X'^)h = 0, where 

B = M'''^di = d^-M'^"". 

Proof. It is clear that the muhiphcation operator M*^^ commutes with the 
Dirac operator ■ Moreover, the operator M'^'^'^ anti-commutes with the 
operator M^. Indeed, as n is even, we have consecutively 

= {-ly-^gfieN 
= -M'^^{Mlg) . 

Furthermore, we have 
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since 

= _(^_l)«("+l)/2 

= 1 , 

where we have used the assumption that — is odd. 
In view of the above, we find that 



= 4 ■ (19) 



Statement (i) follows now easily, since 



where 

A = M'^^dx- 

Equation (fT9|) leads in an analogous way to (ii). □ 
Next we prove the following lemma. 

Lemma 7.1 Let — be odd and f a vector-valued function. Then one has 

(i) ker(A ± A) = ker [d^ - Ml^^^'^^) 
and 

(li) ker (5 ± A) = ker [d^ + Ml^^^^^) . 

Proof. We restrict ourselves to the proof of statement (i), the proof of (ii) 

being similar. 

We have consecutively 

(yl±A)^ = ^ M^"^(a^-M^)^±Ac/ = 

^ {d^- Ml-)g ± M'^^ \g = Q 

^ d^g - g{f_ =F Mcn) = 

^ {d^-Ml'^^'^^)g = , 
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which proves (i) . □ 

We now arrive at the main result of this section. 

Theorem 7.1 Let — he odd and f a vector-valued function. Then the folio- 
wing decompositions hold: 

(i) ker(A2 - A^) = ker (6»^ - Mt^^'"^) ker {d^ - Mir^'"^) . 
and 

(ii) ker - A^) = ker (6»^ + Ml-^'^^) © ker (6>^ + ML+>^'^^) . 
Proof. Again we restrict ourselves to the proof of (i) . As 

{A^ - X')g ^(A-X){A + X)g ^(A + X){A- X)g , 
it is clear that 

ker(A - A) + ker(A + A) C ker(A2 - A^) . 

The converse inclusion 

ker(A2 - A^) c ker(A - A) + ker(A + A) 

also holds. This can be seen as follows. Take g e ker(A^ — A^) . By definition 
we then have 

{A' - X')g ^(A-X){A + X)g = (A + A)(A - A)^ = 
and consequently 

{A + X)g e ker(^ - A) and {A - X)g e ker(A + A) . 
Now, decomposing g as: 

we indeed obtain that g e ker(A — A) + ker(A + A) . 
Moreover, as A 7^ 0, we also have 

ker(A - A) n ker(yl + A) = {0} , 

hence we have shown that 

ker(A2 - A^) = ker(>l - A) ker(^ + A) . 

By means of Lemma 7. 1 this proves the theorem. □ 
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Tl 

Corollary 7.1 Let — he odd and f a vector-valued solution of the Clifford 

Riccati equation. Then any scalar-valued solution (p of the Schrddinger equa- 
tion 

(_A„ - v)(j) = (20) 

can he uniquely represented as cj) = g + h , where g and h are Clifford algehra- 
valued solutions of the equations 

{d^ - M/+^*^'^) g = and {d^ - Ml-^'"^) h = 

respectively. 

Proof. From Proposition 3.1, we know that the Schrodinger equation fl20l) is 
equivalent with 

d[d-f<p = (21) 
where / is a solution of the Clifford Riccati equation. Moreover, taking this 
solution / vector- valued and — odd. Proposition 7.1 implies that fl2T]) can be 
rewritten in the form 

(A" - A')^ = 

hence 

e ker(A2 - A^) = ker {d^ - M/+^*^'^) © ker {d^ - Ml'^'^'') . 
Consequently can be uniquely represented as 

(j) = g + h 

with 

g e ker (5^ - Mt^^'"'') and h G ker - Mir^'^^) . □ 

Tl 

Corollary 7.2 Let — he odd and f a vector-valued solution of the equation 

u = dj-f . 

Then any scalar-valued solution (p of the equation 

(-A„ + m)0 = A20 (22) 

can he uniquely represented as (p = g + h , where g and h are Clifford algehra- 
valued solutions of the equations 

{d^ + M/-^*^'^) g = and {d^^ + M=^-+^''=^) h = 

respectively. 
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Proof. Using equation ([T^ to see that (1221) can be rewritten as 

the proof runs along the same lines as the proof of the previous corollary. □ 
Remark 7.1 If we take 

f = /iCl + /2e2 + . . . + fn-lCn-l 

instead of a vector-valued f and Cn instead of the pseudo- scalar cn , all results 

n 

of this section hold for arbitrary dimension, i.e. without the restriction — odd. 
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